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ABSTRACT 


Potential  vorticity  “thinking”  has  been  used  to  infer  the  balanced  dynamics  of 
predominately  synoptic-scale  weather  features.  This  is  accomplished  through  the  use  of 
the  invertibility  principle  in  which  the  wind  and  mass  fields  can  be  retrieved  from  the 
three-dimensional  distribution  of  potential  vorticity  (PV).  The  most  common  use  of  PV 
thinking  has  been  applied  to  various  forms  of  synoptic-scale  cyclogenesis  and  the 
sensitivity  of  certain  weather  features  had  on  the  subsequent  development. 

Utilizing  the  non-linear  balance  PV  inversion  developed  by  Davis  and  Emanuel 
(1991)  this  study  extends  the  current  use  of  PV  diagnostics  to  atmospheric  features  and 
motions  on  the  order  of  the  meso-  and  storm-scale.  This  study,  in  essence  becomes  a 
feasibility  study  to  examine  whether  a  primarily  synoptic-scale  diagnostic  tool  can  be 
applied  to  much  smaller  scales  and  inherently  more  complex  fields  of  motion  that  occur 
within  a  supercell  thunderstonn.  An  idealized  supercell  simulation  from  the  ARW  will 
be  used  to  examine  the  low  level  thunderstorm  dynamics  from  a  PV  perspective. 

The  results  show  promise  that  the  PV  diagnostic  can  be  applied  to  thunderstorm 
dynamics  given  the  qualitative  results  presented  here.  Refinements  are  necessary  to 
improve  the  quantitative  accuracy  of  this  technique. 
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I.  INTRODUCTION 


Thunderstorm  forecasting  has  and  will  continue  to  be  a  primary  concern  for 
meteorologists  due  to  the  varied  nature  of  convection  and  the  weather  hazards  these 
storms  produce.  Storms  can  produce  a  variety  of  weather  phenomena  that  could  be 
considered  severe  (Johns  and  Doswell  1992)  and  represent  a  significant  hazard  to  the 
public  (Brooks  et  al.  1994a).  The  most  commonly  recognized  threats  being  heavy  rain, 
large  hail,  damaging  winds,  and  tornadoes.  Forecasters  frequently  use  pattern  recognition 
as  a  primary  tool  to  forecast  favorable  environmental  conditions  related  to  the 
development  of  severe  convective  storms.  The  fundamental  components  necessary  for 
severe  convection  are  (i)  sufficient  low  level  moisture,  (ii)  conditional  instability  and  (iii) 
a  trigger,  or  mechanism  to  lift  wann,  moist  low  level  air  to  its  level  of  free  convection 
(LFC)  (Doswell  1987,  Johns  and  Doswell  1992).  A  fourth  condition,  vertical  shear  of  the 
horizontal  winds,  helps  determine  the  convective  mode  after  thunderstorm  initiation  has 
begun  (Gold  and  Nielsen-Gammon  2008a,  part  one  of  four  part  series,  hereafter  GNG1). 

While  questions  remain  about  the  link  between  the  synoptic-scale  environment 
and  the  mode  and  intensity  of  convection,  it  has  been  relatively  agreed  upon  that  the 
synoptic-scale  environment  creates  conditions  favorable  for  severe  weather  (GNG1, 
Johnson  and  Mapes  2001).  Forecasters  usually  expect  an  increase  in  severe  weather 
probabilities  in  the  presence  of  an  upper  level  trough.  However,  there  is  still  continued 
difficulty  in  determining  the  most  common  type  of  severe  weather  that  will  occur  even 
when  forecasters  correctly  anticipate  a  severe  weather  outbreak  (Johns  and  Hart  1993, 
GNG1).  Forecasters  diagnose  environmental  regions  supportive  of  severe  weather  but 
the  importance  of  synoptic-scale  flow  and  subsynoptic-scale  processes  is  still  difficult  to 
detennine  (Doswell  1987,  GNG1).  This  inherent  difficulty  in  forecasting  severe  weather 
arises  from  the  multitude  of  mesoscale  processes  that  work  to  initiate  convection. 

Johnson  and  Mapes  (2001)  describe  the  mesoscale  processes  that  work  to  (i) 
precondition  the  environment,  (ii)  initiate  convection,  and  (iii)  provide  convective 
feedback  to  the  environment.  The  convectively  generated  cold  pool  is  a  mesoscale 

feature  that  can  act  in  each  of  the  three  ways.  It  is  then  logical  to  assume  that  the  cold 
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pool  plays  an  important  role  in  thunderstonn  dynamics.  The  cold  pool  fonns  when 
evaporatively  cooled  air  in  the  downdraft  reaches  the  ground  and  spreads  out  horizontally 
(Engerer  et  al.  2008).  Recent  attention  by  Engerer  et  al.  (2008)  has  been  given  to  the 
study  of  surface  characteristics  of  colds  pools  from  mesoscale  convective  systems  (MCS) 
using  the  Oklahoma  Mesonet  (see  Brock  et  al.  1995  for  a  description).  The  results  found 
by  Engerer  et  al.  (2008)  showed  a  substantial  drop  in  potential  temperature  with  the 
passage  of  the  gust  front  that  represents  the  leading  edge  of  the  cold  pool.  This 
temperature  deficient  could  allow  the  cold  pool  to  be  characterized  as  a  negative  PV 
anomaly. 

The  concept  of  PV  invertibility  has  made  it  possible  to  retrieve  the  wind  and  mass 
fields  from  a  given  distribution  of  PV  based  on  an  appropriately  chosen  balance 
constraint  (e.g.  Hoskins  et  al.  1985  hereafter  HMR85,  Bluestein  1993).  This  technique 
has  been  successfully  applied  to  studies  on  extra-  and  tropical  cyclogenesis  (e.g.  Davis 
and  Emanuel  1991,  hereafter  DE91,  Davis  et  al.  1996,  Huo  et  al.  1999)  and  more 
recently,  has  been  applied  to  the  severe  convective  environment  (e.g.  Roebber  et  al.  1999, 
Gold  and  Nielsen-Gammon  four  part  series  2008)  to  diagnose  mid  tropospheric 
shortwave  troughs.  These  studies  have  used  PV  inversion  to  both  quantify  and 
qualitatively  describe  the  sensitivity  a  certain  upper  level  feature  has  on  a  particular 
weather  event.  Hoskins  et  al.  (1985)  and  Thorpe  (1985)  showed  that  atmospheric  PV 
anomalies  have  a  surface  equivalent,  namely,  potential  temperature  anomalies. 

The  nonlinear  balance  (NLB)  PV  inversion  developed  by  DE91  and  supercell 
simulations  from  the  Weather  Research  and  Forecasting  (WRF)  model  version  three  will 
be  used  in  this  study.  While  most  PV  studies  are  most  commonly  applied  to  synoptic- 
scale  processes,  this  study  hopes  to  extend  the  current  use  of  PV  diagnostics  toward 
mesoscale  processes  since  it  is  those  processes  that  play  a  major  role  in  initiating 
convection.  The  cold  pool  is  an  important  mesoscale  feature  for  the  development  of 
convection  and  simultaneously  be  represented  as  a  surface  temperature  anomaly.  The 
edge  of  the  cold  pool  is  also  important  as  horizontal  vorticity  is  tilted  into  the  vertical  and 
ingested  into  the  thunderstorm’s  updraft.  This  study  becomes,  in  essence,  a  feasibility 
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exercise  to  examine  whether  PV  inversion  can  be  applied  to  mesoscale  and  even  storm- 
scale  processes  with  a  focus  on  low  level  features  that  affect  convective  development. 

The  remaining  sections  of  the  thesis  is  organized  as  follows;  background 
information,  data  and  methodology,  data  and  results  and  finally  conclusions  and  future 
research.  The  Background  (Chapter  II)  contains  infonnation  on  the  (i)  WRF  model  and 
Advanced  Research  WRF  supercell  simulation,  (ii)  potential  vorticity,  and  (iii) 
convective  dynamics  and  characteristics  of  cold  pools.  Chapter  III  describes  the 
supercell  simulation  and  the  potential  vorticity  inversion  technique.  The  data  and  results 
(Chapter  IV)  is  divided  into  three  parts  based  on  the  convective  lifecycle;  (i)  early 
development,  (ii)  stonn  splitting,  and  (iii)  supercell  development.  Chapter  V  provides  the 
conclusions  and  areas  of  future  research 
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II.  BACKGROUND 


A.  THE  ADVANCED  RESEARCH  WRF 

The  WRF  model  is  a  numerical  weather  prediction  model  that  contains  two 
different  systems  designed  to  support  operational  and  research  applications  (Skamarock 
et  al.  2008).  The  operational  subset  is  the  WRF-NMM  (Non-hydrostatic  Mesoscale 
Model)  and  is  used  to  produce  the  standard,  daily  operational  weather  forecasts.  The 
other  subset  is  the  WRF-ARW  (Advanced  Research  WRF)  and  is  used  to  support  a  wide 
range  of  research  applications.  The  ARW  is  designed  to  handle  real-data  and  idealized 
forecasts  and  simulations  due  to  its  large  variety  of  physical  and  dynamical  options.  Use 
of  the  WRF  preprocessing  system  (WPS)  is  required  for  real-data  studies.  The  WPS  is  a 
set  of  programs  that  take  meteorological  and  terrestrial  data  and  transforms  them  into  a 
format  suitable  for  input  into  the  ARW.  An  overview  of  WRF  system  components  is 
given  in  Figure  1 . 


Figure  1.  WRF  system  components.  From  Skamarock  et  al.  (2008). 
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The  ARW  is  composed  of  the  ARW  dynamics  solver  and  the  other  compatible 
components  of  the  WRF  system  (e.g.,  physics  packages,  numeric  and  dynamic  options, 
data  assimilation).  The  latest  version,  3,  was  released  in  April  2008.  All  the  infonnation 
presented  on  the  ARW  has  been  taken  directly  from  the  National  Center  for  Atmospheric 
Research  (NCAR)  Technical  note  NCAR/TN-475+STR  (Skamarock  et  al.  2008)  and  the 
online  WRF  modeling  systems  user  guide  (Wang  et  al.  2008)  located  at 
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/.  The  reader  is  referred  to 
these  sources  for  additional  technical  infonnation  not  presented  here.  Some  additional 
references  are  also  provided  on  specific  features  of  the  ARW  model  where  more 
information  can  be  obtained. 

The  ARW  dynamics  solver  integrates  the  compressible,  non-hydrostatic,  Euler 
equations  in  flux  form  (see  Ooyama  1990).  The  equations  used  a  terrain-following 
hydrostatic-pressure  vertical  coordinate  (Laprise  1992)  commonly  referred  to  as  eta- 
levels. 

The  ARW  is  discretized  over  time  and  space.  Temporal  discretization  is 
performed  through  a  time-split  integration  scheme  described  by  Wicker  and  Skamarock 
(2002).  This  approach  prevents  the  meteorologically  insignificant  (high-frequency) 
modes  from  limiting  the  time  step  of  the  third-order  Runge-Kutta  (RK3)  scheme  used  to 
integrate  the  slower,  meteorologically  significant  modes.  Spatial  discretization  uses  the 
Arakawa  C  grid  staggering.  That  is,  velocities  are  staggered  one -half  grid  length  from 
the  thermodynamic  variables  as  shown  in  Figure  2. 
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Figure  2.  Horizontal  and  vertical  grids  for  the  ARW.  From  Skamarock  et  al.  (2008). 

The  ARW  model  physics  have  many  options  available.  The  five  main  categories 
are  (i)  microphysics,  (ii)  cumulus  parameterizations,  (iii)  planetary  boundary  layer  (PBL), 
(iv)  land-surface  model  (LSM),  and  (v)  radiation.  Each  of  these  five  categories  contains 
many  specific  parameterizations  to  choose  from.  There  are  also  many  interactions 
(Table  1)  between  the  parameterizations  and  the  state  variables  and  their  tendencies. 
Physics  drivers  separate  the  physics  section  from  the  dynamic  solvers.  The  physics 
drivers  are  divided  between  a  pre-physics  preparation  step  and  post-physics  modification 
of  variable  tendencies. 
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Table  1.  Physics  Interactions.  Columns  correspond  to  model  physical  processes: 
radiation  (Rad),  microphysics  (MP),  cumulus  parameterization  (CP),  planetary 
boundary  layer/vertical  diffusion  (PBL),  and  surface  physics  (Sfc).  Rows 
correspond  to  model  variables  where  i  and  o  indicate  whether  a  variable  is  input 
or  output  (updated)  by  a  physical  process.  From  Skamarock  et  al.  (2008) 


Rad 

MP 

CP 

PBL 
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io 
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io 

io 

io 

io 

Tendencies 

Water  Vapor 

io 

io 

io 

io 

Cloud 

io 

io 

0 

io 

Prccip 

io 

io 

0 

Surface  Fluxes 
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Heat  Flux 

i 

0 

Moisture  Flux 

i 

0 

Surface  Stress 

i 

0 

B  POTENTIAL  VORTICITY 

1.  Overview  and  the  Invertibility  Principle 

The  simplest  form  of  PV  is  based  on  the  work  of  Rossby  (1940)  and  Ertel  (1942). 
Describing  the  two  dominant  processes  in  the  vorticity  budget,  vorticity  is  created  by 
vertical  stretching  of  vortex  tubes  and  by  the  horizontal  advection  of  absolute  vorticity. 
The  overall  concept  of  potential  vorticity  asserts  that  there  is  a  potential  for  creating 
vorticity  by  changing  the  latitude  of  the  parcel  and  by  adiabatically  altering  the  distance 
that  separates  the  layer  bounded  by  two  isentropic  surfaces  (HMR85).  Thus,  PV  can  be 
thought  of  as  a  measure  of  the  ratio  of  absolute  vorticity  to  the  effective  depth  of  the  layer 
containing  the  vortex  (Holton  1992).  Rossby  (1940)  and  then  Ertel  (1942)  expressed  this 
relationship  as 

PV  =  -g(J  +  k-Vgxv)^  (1) 

op 
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where  £g  =  k -V  exv ,  is  the  isentropic  vorticity,  k  is  the  vertical  unit  vector  and  V0  is  the 
three  dimensional  gradient  operator  in  x,  y,  0  space  (HMR85).  This  relationship  showed 
that  isentropic  PV,  or  Ertel  PV,  is  conserved  in  an  adiabatic,  frictionless  atmosphere  and 
for  balanced  flow  or  unbalanced  flow  (HMR85,  DE91,  Bluestein  1993,  Morgan  and 
Nielsen-Gammon  1998,  hereafter  MGN98). 

Early  work  using  PV  mainly  focused  on  its  significance  as  an  air-mass  tracer. 
Kleinschmidt,  through  a  series  of  papers  in  the  1950’s,  pushed  the  use  of  PV  and  tried  to 
explain  cyclogenesis  using  upper  tropospheric  PV  anomalies  (HMR85).  This  work  led  to 
the  formulation  known  today  as  the  “invertibility  principle”  for  PV,  namely  that  one  can 
diagnose  the  complete  atmospheric  flow  structure  from  the  spatial  distribution  of  PV 
(Kleinschmidt  1950a,  b,  1951,  HMR85,  Bluestein  1993).  The  ability  to  invert  the  PV 
field  and  retrieve  the  wind  and  mass  fields  allows  PV  to  become  a  useful  diagnostic  tool 
for  the  dynamical  processes  that  occur  in  the  atmosphere. 

The  application  of  the  invertibility  principle  requires  additional  information  in 
order  to  obtain  the  individual  fields  of  absolute  vorticity  and  static  stability.  To  perform 
the  inversion,  the  (i)  reference  state  of  the  atmosphere  must  be  specified,  (ii)  an 
appropriate  balance  condition  must  be  imposed,  and  (iii)  the  inversion  must  be  solved 
globally  with  proper  boundary  conditions  (HMR85).  The  use  of  a  proper  reference  state, 
balance  condition  (e.g.  geostrophy,  gradient  wind  balance,  non-linear  balance)  and 
boundary  conditions  allows  us  to  find  the  exact  combination  of  absolute  vorticity  and 
static  stability  associated  with  the  global  PV  field  (HMR85,  Bluestein  1993). 


2.  Distribution  of  PV  in  the  Atmosphere 

Recall  Rossby’s  form  of  potential  vorticity  in  Eq  (1)  where  PV  is  the  product  of 
the  isentropic  vorticity  and  static  stability.  The  standard  PV  can  be  simplified  to 

PV«-gf^-.  (2) 

op 

Substituting  typical  values  associated  with  mid  latitude  synoptic  scale  flow  into  Eq  (2): 


PV  =  -(10  mV'XIO'V2) 


10  K 
10  kPA 


1  kPa 


1 0  kg  ms  “  in 


(3) 


=  1(T6  m2  V1  K  kg-1  =  1  PVU  (4) 

where  a  PVU  is  a  potential  vorticity  unit.  PV  values  approximately  less  than  1.5  PVU  are 
usually  associated  with  tropospheric  air  while  higher  values  are  associated  with 
stratospheric  air  (Bluestein  1993). 

This  boundary  between  tropospheric  and  stratospheric  PV  can  be  characterized  by 
the  dynamic  tropopause.  The  dynamic  tropopause  is  defined  based  on  a  constant  value  of 
PV  instead  of  the  standard  lapse  rate  definition  (MNG98).  Various  values  of  PV  have 
been  used  to  define  the  tropopause  but  the  most  commonly  used  values  range  from  1-3 
PV  units.  This  range  usually  lies  in  the  transition  region  between  the  upper  troposphere 
and  lower  stratosphere.  Due  to  different  thennal  stratification  and  different 
concentrations  of  certain  chemical  constituents  (e.g.,  ozone  induced  stratospheric 
wanning);  the  stratosphere  tends  to  have  higher  values  of  PV  than  the  upper  troposphere 
(MNG98,  Gerard  1998).  At  high  latitudes,  this  region  of  high  PV  values  is  often  referred 
to  as  the  stratospheric  reservoir.  Equatorward  protrusions  of  relative  PV  maxima  indicate 
regions  where  stratospheric  air  is  descending  into  the  troposphere.  Figure  3  shows  the 
average  vertical  distribution  of  Ertel’s  PV  for  July.  As  expected,  the  tropopause  is 
highest  over  the  tropics  and  slopes  downward  over  the  poles.  The  stratospheric  reservoir 
is  also  easily  recognized  over  the  higher  latitudes  and  poles. 
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latitude 


Figure  3.  Vertical  cross  section  of  zonal  mean  of  Ertel’s  potential  vorticity  in  PVU 
(dashed  lines)  and  potential  temperature  in  K  (solid  lines)  for  July;  from  1979- 
1989  initialized  ECMWF  data  at  seven  levels,  twice  daily.  From  Bluestein 

(1993). 


3.  Atmospheric  PV  Anomalies 

The  invertibility  principle  uses  the  idea  that  the  atmosphere  can  be  described  as 
the  summation  of  all  positive  and  negative  PV  anomalies  and  the  wind  field  is  a 
summation  of  the  individual  wind  fields  attributable  from  each  separate  anomaly 
(Bluestein  1993).  A  correlation  should  then  be  made  that  relates  the  PV  field  to  a  more 
commonly  used  meteorological  field.  Figure  4  shows  the  time  evolution  from  12  UTC 
11  March  2006  to  12  UTC  12  March  2006  using  the  11/12Z  GFS  model  run.  Relative 
maxima  PV  values  correspond  to  cyclonic  curvature  and  cooler  air  associated  with 
troughs  and  closed  lows. 
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Figure  4.  Potential  vorticity  (left)  and  corresponding  500  mb  height  field  in  meters 
(right)  from  the  12  UTC  03/11/2006  GFS  model  run.  Top  row:  1 1/12Z 
initialization.  Middle  row:  12  hour  forecast  valid  at  12/00Z.  Bottom  row:  24 
hour  forecast  valid  at  12/12Z.  Only  positive  PVU  values  greater  than  1  PVU  are 

shown. 
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While  only  positive  PVU  values  are  shown,  the  same  inference  can  be  drawn 
from  relative  minima  PV  values.  PV  minima  correspond  to  anticyclonic  curvature  and 
warmer  air  associated  with  ridges  and  blocking  cyclones.  Thus,  one  can  define  relative 
maxima  (minima)  PV  values  as  positive  (negative)  PV  anomalies  (HMR85,  Bluestein 
1993).  A  note  must  be  made  that  the  positive  (negative)  term  refers  to  the  Northern 
Hemisphere  only.  To  be  precise,  a  positive  (negative)  PV  anomaly  should  be  classified 
as  cyclonic  (anticyclonic)  which  holds  true  for  the  Southern  Hemisphere.  Since  our  case 
study  is  focused  in  the  United  States,  positive  and  negative  terms  will  used. 

We  define  PV  as  the  product  of  absolute  vorticity  and  static  stability  (Eq  1).  A 
positive  PV  anomaly  would  then  be  more  stable  than  a  negative  PV  anomaly  assuming 
absolute  vorticity  being  equal  for  both  anomalies.  Given  this  difference  in  static  stability, 
a  positive  PV  anomaly  would  see  the  isentropes  in  a  vertical  cross  section  bend  toward 
the  center  of  the  anomaly.  This  would  increase  the  static  stability  of  the  anomaly  itself 
while  decrease  the  stability  above  and  below  the  anomaly  due  to  a  weakening  of  the 
potential  temperature  gradient.  Figures  5  and  6  show  the  vertical  cross  sections  of  a 
positive  and  negative  PV  anomaly  respectively. 


Thorpe  (1985)  used  the  non-linear  polar  coordinate  fonn  of  the  PV  equation 


a 

(  1  d(rV )) 

,  C*  1  5 

{  f,x  aG 

dr 

r  dr  j 

cr*  g  d6 

l K(P) 

df  *dP 

+  — =  cr  - 

dr  dr 


(5) 


to  solve  for  idealized,  circularly  symmetric  positive  and  negative  PV  anomalies  where: 


K(P)  =  J- 

op 


C, 


f  \R/CP 

P_ 

Uoy 


,  ,  2V 

floc-f  + - • 


(6) 


(7) 


The  parameter  floc  is  the  absolute  local  rotation  rate  where  /  is  from  the  Earth’s  rotation 


and  2V  /  r  is  the  relative  rotation  around  the  origin.  V  is  the  wind  speed,  £a0  is  the 

absolute  vorticity  on  an  isentropic  surface,  and  a*  is  a  measure  of  the  static  stability 
(Bluestein  1993).  The  full  derivation  of  Eq  (5)  can  be  found  in  Bluestein  (1993).  The 
schematic  results  of  Thorpe’s  (1985)  solutions  are  depicted  in  Figure  7a  for  a  positive  PV 
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anomaly  and  Figure  7b  for  a  negative  PV  anomaly.  For  more  detail  on  the  computational 
method  used,  see  Thorpe  (1985)  and  Bluestein  (1993).  The  observed  structures  of 
atmospheric  anomalies  are  accurately  depicted  in  Thorpe’s  (1985)  solutions  by 
comparing  Figures  5  and  6  with  Figure  7. 

The  positive  PV  anomaly  is  characterized  by  a  relative  high  static  stability  and 
surrounded  by  a  cyclonic  circulation.  The  static  stability  above  the  positive  anomaly  is 
relative  low.  The  negative  PV  anomaly  is  characterized  with  a  relative  low  in  static 
stability  and  an  anticyclonic  circulation.  The  atmosphere  has  high  static  stabilities  above 
and  below  the  negative  anomaly. 

Bluestein  (1993)  and  HMR85  used  Thorpe’s  (1985)  results  to  show  that  vorticity 
and  static  stability  are  positively  correlated  with  PV  anomalies  on  isentropic  surfaces. 
This  is  also  the  true  for  individual  anomalies  of  vorticity  and  static  stability.  Thermal 
wind  balance  suggests  that  a  positive  PV  anomaly  caused  by  vorticity  alone  would  cause 
the  isentropes  to  squeeze  together  toward  the  center  of  the  anomaly  thus  increasing  static 
stability.  The  same  reasoning  also  holds  for  a  PV  anomaly  composed  entirely  of  a  static 
stability  anomaly.  Radial  pressure  on  the  isentropes  would  have  to  correspond  to  an 
increase  in  vertical  shear  and  thus  vorticity  around  the  anomaly  (Bluestein  1993). 
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MON  Mai  17  2008  0000V000  AVN  GFS  1-7  Scalar  Normal  Wind 

MON  Mar  17  20080000V000  AVN_GFS_l-7  Poiencial  Temperature 


Figure  5.  Vertical  cross  section  through  an  upper  level  positive  PV  anomaly  at  00  UTC, 
March  17,  2008.  Cross  section  taken  along  a  north-northwest  to  south-southeast 
line  from  San  Francisco,  Ca  to  Cabo  San  Lucas,  MEX.  Potential  temperature  in 
K  (blue  contours);  scalar  normal  winds  (red  contours).  Tropopause  is  indicated 

by  thick,  black  line. 


TUE  Mar  18  2008  12O0V0O0  AVN  GFS_l-7  Scalar  Normal  Wind 

TUE  Mar  18  2008 120QV000  AVN_GFS_I-7  Potential  Temperature 


Figure  6.  As  in  Figure  5  but  for  a  vertical  cross  section  through  an  upper  level  negative 
PV  anomaly  at  12  UTC,  March  18,  2008.  Cross  section  taken  along  a  northwest 
to  southeast  line  from  North  Platte,  NE  to  Little  Rock,  AR. 
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Figure  7.  The  cross  sectional  view  of  the  schematic  results  from  Thorpe’s  solutions  to 
Eq  (2.5)  for  a  (a)  positive  PV  anomaly  and  (b)  negative  PV  anomaly.  Tropopause 
(thick  solid  line);  center  of  anomaly  located  at  2500  km  near  the  tropopause;  the  + 
and  -  indicate  the  location  of  maximum  normal  wind  flow  into  and  out  of  the 
cross  section  plane,  respectively.  From  Bluestein  (1993) 
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4.  Surface  Potential  Temperature  Anomalies 

Solutions  for  Eq  (5)  were  also  computed  by  Thorpe  (1985)  for  idealized, 
circularly  symmetric  anomalies  in  surface  potential  temperature  with  no  upper  level  PV 
anomaly  present.  Schematic  results  are  shown  in  Figure  8. 

As  with  atmospheric  PV  anomalies,  a  correlation  can  be  seen  between  warm  and 
cold  surface  potential  temperature  anomalies  with  vorticity  and  static  stability.  A  warm 
surface  anomaly  has  higher  static  stability  as  isentropes  bend  downward  toward  the 
anomaly  and  bunch  at  the  surface.  Assuming  an  unperturbed  upper  level  flow,  the  warm 
core  anomaly  must  induce  a  cyclonic  vortex  based  on  anticyclonic  thermal  wind  shear 
present.  This  same  reasoning  can  also  be  applied  to  a  cold  surface  potential  temperature 
anomaly.  However,  the  cold  anomaly  is  not  confined  to  the  surface  but  extends  upward 
like  a  cold  dome.  This  is  due  to  the  weakened  static  stability  as  the  isentropes  bend 
upward  away  from  the  anomaly.  An  anticyclonic  circulation  occurs  around  the  cold 
anomaly  due  to  cyclonic  thennal  wind  shear  (Thorpe  1985,  HMR85,  Bluestein  1993). 

Based  on  the  results  from  Thorpe  (1985),  a  wann  (positive)  potential  temperature 
anomaly  at  the  surface  is  equivalent  to  a  positive  PV  anomaly  and  a  cold  (negative) 
potential  temperature  anomaly  is  equivalent  to  a  negative  PV  anomaly. 
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8.  The  cross  sectional  view  of  the  schematic  results  from  Thorpe’s  solutions  to 
Eq  (2.5)  for  a  (a)  warm  surface  potential  temperature  anomaly  and  (b)  cold 
potential  temperature  anomaly.  Tropopause  (thick  solid  line);  center  of  anomaly 
located  at  2500  km  near  the  ground;  the  +  and  -  indicate  the  location  of  maximum 
nonnal  wind  flow  into  and  out  of  the  cross  section  plane,  respectively.  From 

Bluestein  (1993) 


5. 


PV  Influences  on  the  Severe  Environment 


The  sensitivity  of  atmospheric  variables  can  be  examined  with  the  use  of  PV 
invertibility.  Previous  studies  (e.g.  Dickinson  et  al.  1997,  Huo  et  al.  1999,  Roebber  et  al. 
2002)  have  used  a  concept  called  PV  surgery  to  understand  how  a  specific  PV  anomaly 
affected  a  particular  weather  event.  PV  surgery  involves  removing  a  specific  PV 
anomaly  and  inverting  the  remaining  PV  field  to  quantify  the  effects  of  the  anomaly  in 
question.  In  order  to  use  PV  invertibility  to  diagnose  severe  convective  weather,  the 
convective  parameters  must  be  linked  to  the  balanced  and  unbalanced  wind  and  height 
fields  (GNG1). 


Storm-relative  enviromnental  helicity  (SREH),  vertical  shear  (SHR),  bulk 
Richardson  number  shear  (BRNSHR)  and  convective  available  potential  energy  (CAPE) 
are  some  of  the  most  common  parameters  associated  with  severe  convection.  These 
parameters  can  be  calculated  once  the  wind  and  height  fields  have  been  retrieved  from 
the  PV  inversion  (GNG1)  as  shown  below: 


SREH=  J  (V-C)-VxV dz  =  J  k 

sfc  sfc 


dV 
x - 

dz 


■(V-C)dz 


SRH=  V6-VPBL 


(8) 

(9) 


BRNSHR  =  0.5  x  (u2  +  V 2 ) 


(10) 


lnpo 

CAPE=  J  Rj(T^ 

ln  Plnb 


(11) 


where  z  is  the  assumed  flow  depth,  C  is  the  storm  motion  vector  calculated  using  the 
method  derived  by  Bunkers  et  al.  (2000).  V6  is  the  6  km  wind,  V  is  the  wind  in  the 

planetary  boundary  layer  (PBL),  and  u  and  v  are  the  component  differences  between  the 
mass  weighted  winds  in  the  0-6  km  layer  and  the  PBL.  The  equation  for  CAPE  consists 
of  the  parcel’s  virtual  temperature,  T  ,  the  parcel’s  original  level,  p() ,  and  the  parcel’s 

level  of  neutral  buoyancy,  pLNB  (Brooks  et  al.  1994a,  GNG1).  Operational  forecasters 
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use  these  parameters  to  quantify  the  buoyancy  and  wind  shear  which  are  then  used  to 
determine  the  mode,  intensity,  and  longevity  of  thunderstorms  (Weisman  and  Klemp, 
1982). 

Nielsen-Gammon  and  Gold  (2008,  part  two  of  four  part  series,  hereafter  NGG2) 
demonstrated  the  effects  of  upper  level  PV  changes  on  the  convective  parameters  using  a 
positive  idealized,  circularly  symmetric  PV  anomaly.  As  seen  in  Figure  9,  the  results 
were  consistent  with  those  of  Thorpe  (1985)  with  the  anomaly  decreasing  in  magnitude 
as  the  horizontal  and  vertical  distance  increase.  Applying  these  results  to  the  severe 
convective  environment,  the  downward  bending  of  isentropes  above  the  anomaly  shows  a 
warming  of  that  layer.  Conversely,  the  upward  bending  of  the  isentropes  below  the 
positive  PV  anomaly  corresponded  to  a  cooling  of  the  atmosphere  since  potential 
temperature  increases  with  height.  This  thermal  structure  is  similar  to  the  cold-core 
trough  described  by  Thorpe  (1986)  and  suggests  that  the  region  beneath  and  near  the 
anomaly  has  a  higher  convective  instability  (GNG1). 

Studies  done  by  Jukes  and  Smith  (2000)  and  MNG98  showed  that  CAPE  is 
increased  and  convective  inhibition  is  decreased  directly  under  and  near  the  positive  PV 
anomaly.  This  change  in  the  convective  environment  is  due  to  lower  values  of  potential 
temperature  being  pulled  upward  toward  the  anomaly  causing  a  cooling  in  that  section  of 
the  atmospheric  column.  Also  seen  in  Figure  9  is  an  increase  in  the  cyclonic  circulation 
with  height  below  the  anomaly.  This  increase  in  winds  with  respect  to  height  would 
increase  the  vertical  shear  below  and  near  the  anomaly  and  aid  in  developing  an 
environment  conducive  for  rotating  updrafts  (NGG2). 

The  amplification  of  the  PV  anomaly  would  intensify  the  effects  listed  above  as 

the  temperature  perturbations  increase  above  and  below  the  anomaly.  The  atmospheric 

column  below  and  near  the  intensified  anomaly  would  experience  an  increase  in  the 

upward  bending  of  the  isentropes  thus  increasing  the  cooling.  The  cooling  below  the 

anomaly  would  act  in  the  same  manner  as  cold  air  advection  and  decrease  the  lapse  rates 

in  the  mid  levels  of  the  atmosphere.  This  would  increase  CAPE  and  decrease  CIN  in  the 

column  beneath  and  near  the  anomaly  and  increase  the  probability  of  deep  moist 

convection.  The  cyclonic  circulation  around  the  PV  anomaly  would  also  strengthen  in 
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response  to  the  amplification  of  the  anomaly.  This  would  likely  lead  to  an  increase  of 
vertical  shear  in  the  convective  environment.  Large  scale  assent  would  also  be  affected 
as  air  flows  up  and  along  the  isentropes.  The  upward  bending  of  the  isentropic  surfaces 
would  increase  vertical  motion  beneath  and  downstream  of  the  anomaly.  Rising  air 
would  cool  adiabatically  and  further  cool  the  layer  beneath  the  anomaly  and  add  to  the 
increased  CAPE  and  decreased  CIN  (NGG2). 


^ 


Figure  9.  Idealized  schematic  of  an  isolated  positive  (cyclonic)  spherical  PV  anomaly. 

Background  state  is  uniform  and  the  vertical  axis  is  scaled  by  No/fo.  From  NGG2. 

C.  CONVECTIVE  DYNAMICS  OF  SEVERE  THUNDERSTORMS 
1.  Overview 

The  first  conceptual  model  of  a  thunderstorm  was  derived  from  data  collected 
over  Florida  during  the  Thunderstorm  Project.  This  project  took  place  in  the  post  World 
War  II  years  as  a  joint  venture  between  the  Weather  Bureau,  National  Advisory 


21 


Committee  for  Aeronautics  ,U.S.  Air  Force  and  U.S.  Navy  to  investigate  the  dynamics  of 
thunderstorms  (Byers  and  Braham  1948,  Doswell  2001).  While  the  conceptual  model 
produced  by  Byers  and  Braham  (1948)  has  changed  with  time  and  continuing  research, 
the  concept  of  the  “thunderstorm  cell”  remains  today  as  the  basic  organizational  structure 
for  convection  (Doswell  2001),  namely  the  formation  of  a  single  updraft  and  downdraft 
pair. 

Numerous  studies  have  detailed  the  severe  convective  environment  and  the 
structure  of  supercell  thunderstorms.  The  need  to  accurately  predict  the  formation  of  and 
subsequent  motion  of  supercells  has  provided  50  plus  years  of  research  into  the  topic. 
The  continued  research  on  severe  thunderstorms  has  also  been  guided  by  the  fact  that 
most  supercells  produce  some  type  of  severe  weather  (e.g.  tornadoes,  large  hail, 
damaging  winds  or  flashfloods)  at  some  point  during  the  life  cycle  of  the  storm  (Burgess 
and  Lemon  1991,  Bunkers  et  al.  2000). 

Following  the  terminology  presented  by  Moller  et  al  (1994)  and  Bunkers  et  al 
(2000),  a  supercell  is  defined  as  a  convective  storm  that  contains  a  mesocyclone  or 
antimesocyclone  with  (i)  the  absolute  value  of  vertical  vorticity  greater  than  or  equal  to 
10'  s'  ,  (ii)  persist  at  a  minimum  on  the  order  of  tens  of  minutes,  and  (iii)  persist  through 
at  least  one  third  of  the  storms  depth.  Right  and  left  movers  are  classified  by  the 
direction  they  move  relative  to  the  vertical  shear.  Right  (left)  movers  are  mesocyclones 
(antimesocyclones)  that  travel  to  the  right  (left)  of  the  vertical  wind  shear. 

The  dynamics  of  the  severe  thunderstonn  environment:  how  thunderstorms 
develop  and  progress  through  their  convective  life  cycles,  has  been  well  documented 
(e.g.,  Lemon  and  Doswell  1979,  Weisman  and  Klemp  1982,  Rotunno  and  Klemp  1985, 
Moller  et  al.  1994)  along  with  the  characteristic  differences  between  tornadic  and 
nontomadic  storms  (e.g.,  Brooks  et  al.  1994a,  1994b).  Since  the  goal  of  this  paper  is  to 
analyze  the  formation  of  rotating  updrafts  from  a  PV  perspective,  the  focus  should  be 
placed  on  those  environmental  conditions  and  dynamics  that  create  mesocyclones. 
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2. 


The  Role  of  the  Cold  Pool  in  Convection 


As  already  noted  previously,  the  cold  pool  can  affect  the  convective  environment 
in  multiple  ways  (Johnson  and  Mapes  2001).  Byers  and  Braham  (1948)  and  Newton  and 
Newton  (1959)  noted  that  the  cold  pool  was  instrumental  in  the  formation  of  new 
convection  along  the  gust  front  or  outflow  boundary.  The  airmass  discontinuity 
represented  by  the  gust  front  acts  in  much  the  same  way  as  a  cold  front.  In  some  cases, 
conditionally  unstable  environmental  air  is  lifted  to  its  LFC  along  the  gust  front  to  form 
new  convective  cells.  This  is  the  primary  process  for  the  sustenance  of  multicell 
thunderstorms  (Weisman  and  Klemp  1986,  Engerer  et  al.  2008).  Johnson  and  Mapes 
(2001)  also  describe  the  way  cold  pools  precondition  the  environment.  In  these  cases 
forced  lifting  of  conditionally  unstable  air  is  unable  to  initiate  convection  by  itself.  The 
development  of  new  convection  then  requires  additional  forcing.  This  can  occur  as  the 
cold  pool  intersects  other  low  level  boundaries  (e.g.,  other  outflows,  terrain,  larger-scale 
fronts)  to  produce  sufficient  lift  for  convection.  Purdom  (1976)  showed  using  satellite 
data  from  the  Geostationary  Operational  Environment  Satellite  (GOES)  that  intersecting 
boundaries  are  often  a  focus  for  the  intensification  of  convection.  On  larger  scales, 
multiple  cold  pools  can  merge  to  cover  hundreds  of  square  kilometers  and  help  to 
generate  new  convection  hundreds  of  miles  from  the  origin  point  of  the  cold  pool 
(Engerer  et  al.  2008).  The  larger  and  stronger  the  difference  between  the  cold  pool  and 
nearby  environment,  the  longer  such  contrast  can  exist  and  spark  new  convection  hours 
after  the  initial  convection  has  ended  (Doswell  2001,  Johnson  and  Mapes  2001). 

The  cold  pool  is  also  an  important  feature  for  the  generation  of  deep,  moist 
convection  (Rotunno  et  al.  1988)  especially  when  combined  with  vertical  shear  of  the 
environmental  winds  (Szeto  and  Cho  1994).  Not  only  does  the  outflow  boundary  act  to 
trigger  new  convective  cells  to  sustain  the  propagation  of  squall  lines  and  other  types  of 
MCSs,  but  increased  vertical  shear  and  convergence  along  the  gust  increases  chances  of 
severe  weather  (Maddox  et  al.  1980).  Studies  by  Maddox  et  al.  (1980)  and  Markowski  et 
al.  (1998)  showed  that  pre-existing  cold  pools  can  lead  to  the  development  of  strong 
tornadoes  with  the  passage  of  supercells  over  the  cold  pool  (Engerer  et  al.  2008).  On  the 
smaller,  convective  scales,  cold  pools  generated  by  the  forward-flank  and  rear-flank 

23 


downdrafts  (FFD  and  RFD  respectively)  have  been  thought  to  play  an  important  part  in 
tornadogenesis.  Brooks  et  al.  (1994b)  state  that  the  occurrence  of  tornados  in  supercells 
is  closely  linked  to  the  formation  of  a  low  level  mesocyclone  created  from  baroclinically 
generated  vorticity  in  the  evaporatively  cooled  downdrafts  (Rotunno  and  Klemp  1985, 
Davies- Jones  and  Brooks  1993). 

While  cold  pools  have  been  intimately  linked  to  the  convective  life  cycle,  a 
balancing  act  must  occur  between  the  cold  pool  and  surrounding  environment  to  promote 
continued  convection.  Thorpe  et  al.  (1982)  proposed  that  favorable  conditions  for  strong, 
long-lived  convection  occurred  when  low  level  shear  of  the  environment  was  sufficiently 
strong  to  prevent  the  cold  pool  from  moving  away  from  the  cell  (Rotunno  et  al.  1988). 
Otherwise,  the  cold  pool  could  propagate  away  from  the  cell  and  cut  off  the  warm,  moist 
environmental  air  feeding  the  convective  updraft  and  causing  the  thunderstorm  to  decay 
(Engerer  et  al.  2008). 


3.  Characteristics  of  Cold  Pools 


The  cold  pool  originates  inside  the  downdraft  as  evaporatively  cooled  air  reaches 
the  ground  and  spreads  out  along  the  surface.  The  accumulation  of  this  cooler  air  creates 
a  hydrostatic  increase  in  surface  pressure  that  is  a  function  of  temperature  and  cold  pool 
depth  (Wakimoto  1982,  Engerer  2008).  This  area  of  higher  pressure  is  commonly 
referred  to  as  a  mesohigh  (Engerer  2008). 


Bryan  et  al.  (2005)  measured  the  intensity  of  MCS  generated  cold  pools  in  the 
stratiform  precipitation  regions  from  soundings  taken  during  the  Bow  Echo  and  MCV 
Experiment  (BAMEX).  Bryan  et  al.  (2005)  measured  cold  pool  intensity  by  C  derived 
from  the  deep  anelastic  equations: 


C2=-- - -f  (pB)dz, 

p(z  =  0)  J0 


(12) 


where  z  is  the  height,  H  is  the  cold  pool  depth,  p  is  the  density,  and  B  is  the  buoyancy 
(neglecting  hydrometeors): 


B=g 


(  a a  \ 

—=  f-0.61(<7v  -gv) 

V  V  J 


(13) 
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In  Eq  (13),  g  is  the  gravitational  acceleration,  0  is  the  potential  temperature,  and  qv  is 
the  water  vapor  mixing  ratio.  The  approximate  height  of  the  cold  pool  was  calculated  to 
be  the  height  above  ground  where  the  buoyancy  first  reached  0 ms  2.  Over-bars  indicate 
variables  obtained  from  the  environment  conditions  with  other  variables  calculated  from 
inside  the  cold  pool. 

The  cold  pools  studied  during  BAMEX  by  Bryan  et  al.  (2005)  showed  that  most 
of  the  cold  pools  were  deep  and  intense,  with  heights  greater  than  three  kilometers  and 
values  of  C  greater  than  20ms  2 .  The  authors  noted  that  results  of  cold  pool  strength 
and  intensity  from  the  BAMEX  study  are  different  than  values  obtained  from  numerical 
studies.  Bryan  et  al.  (2005)  argue  that  the  results  are  an  accurate  representation  of  the 
environments  during  BAMEX  and  that  the  discrepancies  likely  resulted  from  these 
environmental  differences. 

Engerer  et  al.  (2008)  used  the  Oklahoma  Mesonet  to  identify  surface 
characteristics  the  MCS  cold  pools  for  four  MCS  life-cycle  stages,  (i)  first  storms,  (ii) 
MCS  initiation,  (iii)  mature  MCS,  and  (iv)  MCS  dissipation.  Engerer  et  al.  (2008)  found 
that  average  surface  pressure  rises  from  all  stages  was  around  3  hPa.  They  also  noted 
that  the  average  potential  temperature  deficits  inside  the  cold  pool  ranged  from  9.8  to  5.6 
K  colder  with  the  strongest  deficit  occurring  during  the  first  storms  stage  and  decreasing 
an  average  of  1  K  through  each  successive  stage.  Engerer  et  al.  (2008)  did  observe  a 
mean  potential  temperature  deficit  of  nearly  1 1  K  during  the  daytime  and  early  evening 
observations.  Using  the  low  level  potential  temperature  deficit  associated  with  cold  pools 
and  the  baroclinic  generation  of  vorticity  along  the  gust  front,  this  study  hopes  to  exploit 
the  application  of  PV  inversion  to  examine  the  feasibility  of  using  this  dynamical 
technique  on  smaller-scale  features. 

4.  Thunderstorm  Dynamics 

The  ability  to  detennine  convective  storm  types  has  typically  relied  on  the  use  of 
buoyancy,  usually  measured  by  CAPE,  and  the  vertical  wind  shear  (Moller  et  al.  1994). 
These  two  parameters  have  been  shown  to  be  useful  in  discriminating  between  supercell 
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and  non-supercell  storm  environments  but  not  necessarily  between  tomadic  and  nontomadic 
storms  (Brooks  et  al  1993,  Brooks  et  al.  1994b).  Droegemeier  et  al.  (1993),  with  the  use  of 
numerically  generated  storms,  showed  that  the  storm-relative  environmental  helicity  (Eq  8) 
was  an  indicator  of  the  correlation  between  vertical  velocity  and  vorticity.  Higher  helicity 
environments  are  then  assumed  to  have  a  greater  likelihood  of  supporting  rotating  updrafts. 

A  numerical  modeling  study  by  Weisman  and  Klemp  (1982)  showed  that  low  level 
vertical  shear  in  an  unstable  (positive  CAPE)  environment  was  the  origin  of  midlevel 
rotation.  This  conclusion  was  also  supported  by  theoretical  studies  done  by  Davies-Jones 
(1984)  and  Rotunno  and  Klemp  (1985).  Davies-Jones  (1984)  explained  rotating  updraft 
formation  using  streamwise  vorticity.  Convection  deforms  isentropic  surfaces,  which  causes 
a  tilting  of  the  vortex  tubes  that  lie  along  the  isentropic  surface.  Rotunno  and  Klemp  (1985) 
showed  that  the  origin  of  mid-level  updraft  rotation  was  from  the  horizontal  vorticity  from 
environmental  shear.  The  reader  is  referred  to  these  two  papers  for  more  detailed 
information.  Brooks  et  al.  (1994b)  state,  however,  that  the  occurrence  of  tornados  in 
supercells  is  closely  linked  to  the  formation  of  a  low  level  mesocyclone  created  from 
baroclinically  generated  vorticity  in  the  evaporatively  cooled  downdrafts  (Rotunno  and 
Klemp  1985,  Davies-Jones  and  Brooks  1993). 

While  different  mechanisms  are  responsible  for  the  rotation  found  to  occur  in 
thunderstorms,  they  all  have  a  common  trait.  These  mechanisms,  to  a  high  degree,  are  not 
controlled  by  the  large-scale  environment  but  by  smaller-scale  processes.  Using  the 
definitions  of  large-scale  and  mesoscale  provided  by  Doswell  (1987),  large-scale  features  are 
those  that  can  be  described  by  quasi-geostrophic  processes  and  mesoscale  features  are  those 
that  can  be  described  only  by  considering  the  interactions  of  both  large-scale  and  microscale 
processes.  The  evolution  of  synoptic  scale  features  often  creates  an  environment  favorable 
for  severe  weather  (e.g.  low  level  moisture  and  conditional  instability).  However,  large-scale 
lifting  based  on  the  quasi-geostrophic  omega  equation  is  insufficient  to  lift  a  parcel  in  its  LFC 
on  time  scales  necessary  to  initiate  deep,  moist  convection  thus  it  is  the  mesoscale 
interactions  that  affect  the  occurrence  of  convection  (Doswell  1987).  Consequently,  the  goal 
of  this  paper  is  to  examine  low  level  thunderstorm  dynamics  from  a  PV  perspective  in  an 
effort  to  test  whether  PV  inversion  can  be  applied  to  smaller,  more  complex  atmospheric 
motions. 
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III.  DATA  AND  METHODOLOGY 


A.  ADVANCED  RESEARCH  WRF  MODEL 

1.  Data,  Model  Setup,  and  Methodology 

Data  for  this  analysis  will  be  obtained  using  the  idealized  supercell  simulation 
from  version  three  of  the  ARW.  The  domain  will  use  the  Cartesian  coordinate  system 
with  a  grid  box  160  km  x  160  km  in  both  horizontal  directions  and  20  km  in  the  vertical. 
Horizontal  grid  spacing  is  two  kilometers  and  vertical  grid  spacing  is  500  meters.  The 
simulation  will  be  run  through  three  hours  with  a  12  second  time  step  and  WRF  output 
data  written  at  every  five  minute  forecast  time  intervals. 

The  supercell  will  utilize  the  Purdue  Lin  microphysics  scheme.  This  scheme 
contains  six  hydrometeor  classes:  water  vapor,  rain,  snow,  cloud  ice,  cloud  water,  and 
graupel  or  hail.  No  cumulus  scheme  is  used  with  the  two  kilometer  grid  spacing  (Dudhia 
2006)  allowing  for  the  explicit  cumulus  development  with  resolved  updrafts.  The  default 
settings  for  the  supercell  simulation  neglect  both  the  longwave  and  shortwave  radiation 
schemes  and  the  surface-  and  boundary-layer  schemes.  These  default  settings  will  be 
kept  in  order  to  simplify  the  simulation.  The  Coriolis  terms  in  the  model  are  turned  off 
due  to  the  short  (three  hour)  time  period  of  the  simulation. 

The  supercell  simulation  will  be  run  under  ideal  conditions  using  the  default  input 
sounding.  The  supercell  simulation  requires  an  input  sounding  that  contains  any  set  of 
levels  up  to  the  model  top.  The  first  line  of  the  input  sounding  contains  surface  data,  the 
surface  pressure  (hPa),  potential  temperature  (K),  and  moisture  mixing  ratio  (g/kg).  Data 
for  the  input  sounding  contains  the  height  above  the  surface  (in  meters),  the  potential 
temperature  (K),  vapor  mixing  ratio  (g/kg)  and  the  u  and  v  wind  components.  The  model 
then  interprets  and  extrapolates  (if  necessary)  the  data  before  starting  the  simulation.  The 
environmental  winds  of  the  input  sounding  produce  a  quarter  circle  hodograph.  Both  left 
and  right  moving  supercells  are  produced  with  this  input  sounding.  Figure  10  shows  the 
sounding  used  to  initialize  the  supercell  simulation. 
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Figure  10.  Input  sounding  used  to  initialize  the  supercell  simulation 

Using  the  concepts  developed  by  Thorpe  (1985)  and  HMR85  and  focusing  on  the 
low  level  thunderstorm  structure  and  dynamics,  PV  inversion  will  be  used  to  see  if  the 
inverted  wind  and  mass  fields  resemble  the  observed  characteristics  of  thunderstorms. 

2.  PV  Inversion  Technique 

The  ability  to  invert  the  PV  field  requires  the  use  of  a  balance  system  to  relate  the 
wind  and  mass  fields  to  the  distribution  of  PV.  The  balance  system  contains  a  balance 
equation  that  relates  the  streamfunction  to  the  geopotential,  a  vorticity  equation, 
thermodynamic  equation  and  mass  conservation  equation  (MNG98).  Since  the 
convective  environment  is  highly  ageostrophic  in  nature,  a  high  order  balance  constraint 
must  be  used  in  order  to  accurately  retrieve  the  wind  and  mass  fields  from  the  distribution 
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of  PY.  Specifically,  DE91  developed  a  method  for  PV  inversion  using  the  nonlinear 
balance  condition  devised  by  Charney  (1955)  for  Ertel  PV. 

While  PV  diagnostics  using  the  nonlinear  balance  system  have  focused  mainly  on 
synoptic  scale  cyclogenesis,  the  NLB  inversion  should  be  applicable  over  a  large  range  of 
atmospheric  conditions  (Nielsen-Gammon  and  Gold  2008a,  hereafter  NGGa).  The  PV 
framework  provided  by  DE91  makes  use  of  the  NLB  system  in  which  the  irrotational 
component  of  the  horizontal  wind  field  is  small  and  can  be  neglected.  This  creates  a 
small  difference  between  the  NLB  PV  and  the  hydrostatic  Ertel  PV  (NGGa).  For  a  given 
distribution  of  Ertel’s  PV,  the  NLB  PV  inversion  of  DE91  returns  the  following  balanced 
variables;  geopotential  height  and  nondivergent  streamfunction  (DE91,  GNG1).  The 
reader  is  encouraged  to  consult  DE91  for  more  information  on  the  inversion  technique. 


The  PV  inversion  technique  used  here  follows  that  of  DE91.  The  NLB 
(T*,®)  and  potential  vorticity  (q)  equations  are  defined  as 
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with  n  as  the  Exner  function.  The  equations  are  both  scaled  by  appropriate  factors  and 
the  following  equations  used  to  formulate  the  PV  inversion  system  become: 


V2®_V/V'P+v2vF+^ 


/  / 


/ 


a2vp  a2xp 

( a2vp^ 

2\ 

^a*2  dy2 

ydxdy ) 

y 

(14) 


p-q 

gkU 


=  (/  +  V2'P) 


a2® 

arF 


a2vp  a2®  a2vp  a2® 
an  ox  an  ax  anav  anav 


(15) 


The  mass  and  wind  fields  are  linked  to  the  spatial  distribution  of  PV  by  subtracting  Eq 
(14)  and  Eq  (15)  and  by  adding  Eq  (14)  and  Eq  (15).  The  two  equations  used  to  solve  the 
PV  inversion  process  are 
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where  the  geopotential  height  and  streamfunction  fields  are  solved  for  based  on  the 
distribution  of  PV.  The  solution  uses  a  successive  overrelaxation  method  following 
DE91.  The  PV  field  used  in  the  calculations  is  derived  from  the  model  forecasted  wind 
field  and  thermodynamic  structure.  PV  inversion  is  then  used  to  calculate  the  following 
variables:  (i)  u  and  v  wind  components,  (ii)  potential  temperature,  (iii)  geopotential 
height,  and  (iv)  streamfunction  and  streamlines. 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.  EARLY  INITIATION 

The  thermal  bubble  used  to  initiate  convection  is  located  at  the  center  of  the 
domain.  It  is  a  two  degree  Kelvin  temperature  difference  with  a  horizontal  radius  of  10 
kilometers  and  1500  meter  vertical  extent  that  forces  a  parcel  to  its  LFC.  Model  derived 
radar  echoes  are  first  recorded  around  15  minutes  into  the  simulation  associated  with  a 
towering  cumulus  cloud.  Mean  seal  level  pressure  and  950  mb  winds  are  plotted  in 
Figure  1 1  along  with  streamlines. 


Figure  1 1 .  Plot  of  mean  sea  level  pressure  (purple)  and  the  PV  calculated  950  mb  winds 
in  (a)  and  950  mb  streamlines  in  (b)  at  15  minutes  into  the  simulation. 


Developing  convection  is  represented  by  the  low  pressure  in  the  center  of  the 
domain  that  corresponds  to  the  location  of  the  thermal  perturbation.  As  the  rising  column 
of  air  cools  and  condenses  into  a  towering  cumulus  cloud  (Figure  12),  the  PV  derived 
wind  field  begins  to  show  a  low  level  cyclonic  perturbation  around  the  forming  updraft. 
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Figure  12.  Simulated  thunderstorm  in  the  towering  cumulus  stage  20  minutes  into  the 
simulation.  Top  frame  is  looking  north  at  the  vertical  extent  of  the  cloud.  Bottom 
right  is  the  top-down  view  of  the  thunderstorm  and  bottom  left  in  the  derived 
reflectivity.  The  cloud  (in  orange)  is  isosurfaced  at  95%  relative  humidity. 


Recall  the  input  sounding  (Fig.  10)  used  to  initialize  the  simulation  has  westerly 
winds  that  increase  with  height.  The  drier  westerly  winds  are  ingested  into  the  storm 
since  storm  propagation  is  slower  than  the  mid  and  upper  level  winds.  The  drier  winds 
act  to  evaporate  the  rain  and  cloud  water  in  the  storm.  The  evaporatively  cooled  air 
becomes  more  dense  and  falls  to  the  surface  and  spreads  out  horizontally  to  form  the  cold 
pool.  As  the  cold  pool  forms  under  the  storm,  low  level  air  is  forced  upward  into  the 
updraft.  The  baroclinic  region  along  the  edge  of  the  gust  front  is  also  a  region  where 
horizontal  vorticity  is  tilted  in  to  the  vertical.  This  feature  can  be  seen  in  a  cross  section 
through  the  stonn  along  with  the  initial  stonn  rotation  starting  to  develop  in  Figure  13. 
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Figure  13.  Cross  section  of  potential  vorticity  from  975  mb  to  250  mb  through  the 

thunderstorm  at  15  minutes  into  the  simulation.  Two  areas  of  positive  values  are 
noted,  one  at  the  surface  and  one  near  850  mb. 


The  baroclinic  generation  of  vertical  vorticity  can  be  easily  seen  near  the  surface 
in  the  PV  contours  as  well  as  the  beginning  stages  of  storm  rotation.  This  is  seen  by  the 
positive  PV  values  near  850  mb  and  500  mb.  Also  noted  in  the  cross  section  is  the  region 
denoted  by  the  negative  values  in  PV.  As  horizontal  vortex  lines  are  displaced  into  the 
vertical,  a  cyclonic  rotation  is  seen  associated  with  the  updraft  while  an  anticyclonic 
rotation  is  associated  with  downdraft.  The  calculated  PV  field  corresponds  to  these 
features. 

B.  CELL  SPLITTING 

The  original  storm  begins  to  undergo  splitting  between  30  and  45  minutes  into  the 
simulation.  This  split  produces  a  left  and  right  moving  supercell  seen  in  Figure  14.  Right 
movers  produce  deep  cyclonic  circulations  while  left  movers  develop  anticyclonic 
circulations  due  to  increasing  vorticity  that  is  advected  into  the  updraft  (Bluestein  1993). 
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Since  the  left  and  right  moving  storms  have  different  rotation,  an  appropriate  response 
should  be  seen  in  the  PV  distribution  and  the  PV  derived  fields  in  the  regions  around  the 
two  storms  (Figure  15). 


Figure  14.  Model  derived  radar  reflectivities  showing  the  splitting  supercells.  Images 
taken  every  20  minutes.  Radar  sections:  Top  left  at  40  minutes,  top  center  at  60 
minutes,  bottom  left  at  80  minutes  and  bottom  center  at  100  minutes.  The  right 
panels  show  the  thunderstorm  with  relative  humidity  isosurfaced  at  95%  with  the 
top  right  looking  north  and  the  bottom  right  looking  top  down  at  50  minutes  into 

the  simulation. 


The  convection  at  one  hour  is  located  near  the  leading  edge  of  the  cold  pool.  The 

low  level  easterly  winds  have  prevented  the  cold  pool  from  running  ahead  of  the 

thunderstonn  and  cutting  off  the  supply  of  warm  air  from  the  updraft.  This  allowed  for 

the  sustainment  of  each  of  the  supercells.  Figure  15  shows  the  streamline  analysis  at  one 

hour  into  the  simulation.  The  original  storm  now  identified  as  the  left  moving  storm  is 

still  the  dominant  stonn.  The  left  moving  cell  has  taken  on  anticyclonic  flow  as  the  right 

moving  cell  begins  to  develop  along  the  edge  of  the  cold  pool. 

34 


Figure  15.  PV  calculated  fields  at  one  hour  into  simulation  and  at  950  mb  for  (a)  mean 
sea  level  pressure  in  red  and  potential  temperature  in  purple  and  in  (b)  streamlines 

in  purple. 


Cross  sections  taken  though  both  the  left  and  right  moving  storms  further  support 
the  underlying  dynamics.  Figure  16  shows  both  cells  characterized  by  a  nearly  vertical 
shaft  of  PV  associated  with  the  updrafts  in  each  of  the  storms.  The  anticyclonic  updraft 
of  the  left  mover  is  given  by  the  negative  values  in  PV  while  the  cyclonic  updraft  in  the 
right  mover  is  represented  by  positive  PV  values.  The  lower  levels  beneath  both 
thunderstorms  are  characterized  by  strong  positive  values  of  PV  that  signifies  the 
baroclinic  generation  of  vorticity  up  and  over  the  cold  pool  into  the  storms’  updrafts. 
Another  feature  to  note  in  the  right  mover  cross  section  is  the  area  of  negative  PV  values 
located  in  the  lower  levels  on  the  western  flank  of  the  storm.  These  negative  PV  values 
are  most  likely  a  result  of  the  cold  pool  left  in  the  wake  of  the  thunderstonn.  The 
isentropes  are  bulging  upward  in  the  same  manner  as  those  in  Thorpe  (1985)  for  a 
negative  surface  potential  anomaly. 


35 


Figure  16.  Cross  sections  take  through  the  right  mover  (top)  and  left  mover  (bottom). 

The  right  mover  is  characterized  by  upward  bulging  isentropes  (red)  in  the  low 
level  associated  with  the  pool  as  well  as  negative  PV  values  on  the  western 
section  of  the  cold  pool.  Positive  PV  values  correspond  to  a  cyclonic  updraft  and 
large  values  of  PV  along  the  leading  edge  correspond  to  the  generation  of  vertical 
vorticity.  The  left  mover  has  the  same  features  (isentropes  not  shown)  with  the 
exception  of  negative  PV  with  the  anticyclonic  storm  rotation. 
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The  balanced  flow  based  on  the  NLB  PV  inversion  technique  through  the  early 
stages  of  convective  development  and  the  beginning  of  stonn  splitting  have  so  far  been 
qualitatively  accurate  with  other  observational  and  other  numerical  studies  of 
thunderstorms.  These  results  seem  promising  at  least  for  the  beginning  stages.  It  remains 
to  be  seen  however,  how  well  the  NLB  PV  inversion  will  hold  up  as  these  cells  intensify 
and  the  observed  flow  becomes  more  unbalanced. 


C.  SUPERCELL  DEVELOPMENT 

The  right  moving  stonn  continues  to  evolve  as  it  maintains  a  low  level  baroclinic 
vorticity  contribution  to  its  low  level  rotation.  The  right  mover  begins  to  evolve  into 
more  of  the  classic  supercell  as  seen  on  radar.  By  90  minutes  into  the  simulation,  the 
right  mover  has  begun  to  develop  a  “v-notch”  and  the  formation  of  a  pendant  or  hook 
echo  signature  in  the  model  derived  radar  reflectivities  which  suggest  an  intensification 
of  the  thunderstorm.  Figure  17  shows  the  simulated  supercell  at  two  hours  into  the 
simulation. 


Figure  17.  As  in  Figure  12,  but  at  two  hours  into  the  supercell  simulation 
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The  intensifying  supercell  presents  some  problems  for  the  PV  inversion  since  we 
are  using  the  NLB  constraint.  Atmospheric  motions  associated  with  thunderstorms  are 
highly  ageostrophic  so  the  application  of  a  balance  constraint  to  retrieve  an  accurate 
estimate  of  the  storm’s  wind  field  from  the  PV  inversion  might  prove  problematic. 

Figure  18  shows  the  mean  sea  level  pressure  and  wind  field  at  two  hours  into  the 
simulation.  There  are  noticeable  discrepancies  associated  with  the  balanced  wind  speeds 
in  the  vicinity  of  the  two  supercells.  While  the  wind  speeds  are  over-forecasted  near  the 
thunderstorms,  especially  the  right  mover,  the  flow  pattern  around  each  of  the  storms  in 
correct.  The  right  mover  has  cyclonic  flow  and  the  left  mover  has  anticyclonic  flow. 
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Figure  18.  The  balanced  winds  retrieved  from  the  NLB  PV  inversion  at  two  hours  into 

the  simulation  with  wind  barbs  at  950  mb. 


Over-looking  the  gross  errors  from  the  balanced  winds,  the  inversion  process  does 
identify  the  generation  of  new  convection/splitting  cell  formation  associated  with  the 
right  moving  supercell.  Focusing  on  the  right  mover,  the  balanced  streamlines  (not 
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shown  here)  identify  two  cyclonic  rotations  associated  with  two  low  pressures  and  an 
anticyclonic  rotation  separating  the  updrafts. 

To  further  understand  the  possible  reason  behind  the  over- forecasted  wind  speeds, 
Figure  19  shows  the  absolute  vorticity  and  balanced  streamfunction.  High  values  of 
vorticity  are  seen  around  both  thunderstorms.  Focusing  on  the  right  moving  supercell, 
there  is  a  large  value  of  low  level  vorticity  associated  with  the  low  level  cyclonic 
circulation.  This  value  of  absolute  vorticity  is  orders  of  magnitude  larger  than  normal  for 
larger  scale  systems. 


Figure  19.  The  absolute  vorticity  (right)  and  streamfunction  (left)  along  with  streamlines 

(purple)  at  950  mb. 


Recalling  the  Rossby/Ertel  form  of  PV  in  Eq  (1),  gravity  and  the  Coriolis  are 
constants  and  the  static  stability  around  the  right  mover  is  not  changing  enough  to 
generate  this  kind  of  change  in  the  PV  field.  This  leaves  the  relative  vorticity  as  the 
likely  source  of  the  increase  in  PV.  Figure  20  shows  the  three-dimensional  evolution  of 
the  right  mover  beginning  at  80  minutes  into  the  simulation  through  two  hours. 
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Figure  20.  Time  series  of  simulated  tornadic  development.  Moving  from  left  to  right, 
series  starts  at  70  minuets  in  the  simulation  and  ends  at  120  minutes  in  the  bottom 
right.  Clouds  (white)  isosurfaced  at  95%  relative  humidity  with  model  derived 

reflectivity  values  on  the  bottom. 


As  the  right  mover  evolves  and  intensifies  along  the  leading  edge  of  the  cold  pool, 
the  low  level  rotation  is  constantly  being  fed  from  the  low  level  baroclinic  generation  of 
vorticity.  As  the  circulation  becomes  tighter  and  stronger  during  tomadogenesis,  the 
derived  balance  wind  field  must  see  a  response.  The  direct  consequence  to  the  dramatic 
increase  in  relative  vorticity  is  the  unproportional  balanced  wind  field  during  this  period 
of  the  simulation. 

Overall,  the  qualitative  ability  of  the  NLB  PV  inversion  to  handle  storm-scale 
dynamics  was  impressive.  During  the  early  parts  of  convective  development  the 
retrieved  wind  and  mass  fields  closely  resembled  those  observed  in  the  simulation.  The 
later  stages  of  convection  produced  more  of  a  problem  since  the  highly  ageostrophic 
nature  of  convection  puts  strain  on  most  any  balanced  flow  constraint. 
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V.  CONCLUSIONS  AND  FUTURE  RESEARCH 


A.  SUMMARY  AND  CONCLUSIONS 

The  goal  of  this  study  was  to  see  if  PV  inversion  that  utilizes  the  NLB  constraint 
could  be  successfully  applied  to  meso-  and  even  storm-scale  dynamics.  The  results 
presented  here  show  that  there  is  a  good  qualitative  agreement  from  the  balanced  flow  to 
what  is  observed  in  the  thunders tonns.  This  leads  us  to  the  ability  to  characterize  storm- 
scale  features  in  tenn  of  PV  anomalies  to  diagnose  thunderstonn  dynamics. 

Cyclonically  (anticyclonically)  rotating  updrafts  (downdrafts)  or  right  (left) 
moving  stonns  can  be  described  in  terms  of  a  positive  (negative)  PV  anomaly.  The 
surface  cold  pool  which  also  plays  an  important  role  in  convective  development  is  seen  to 
possess  both  aspects  of  positive  and  negative  PV  anomalies.  The  leading  edge  of  the  cold 
pool  is  a  region  of  baroclinic  vorticity  generation  and  thus  a  positive  PV  anomaly,  while 
the  trailing  end  resembles  the  cold  dome  described  by  Thorpe  (1985)  as  a  negative 
surface  potential  temperature  anomaly.  Even  though  the  balanced  flow  was  shown  to  be 
in  qualitative  agreement  with  the  observed  flow,  more  diagnostics  are  needed  to  verify 
the  accuracy  of  this  technique  as  it  has  been  applied  here. 

B.  FUTURE  RESEARCH 

The  choice  of  the  NLB  constraint  seems  to  be  the  biggest  cause  of  error  in  the 
inversion  technique.  While  the  early  time  periods  appear  to  show  some  accuracy,  the 
later  periods  of  the  simulation  have  grossly  over-forecasted  winds  due  to  the  dramatic 
increase  in  low  level  relative  vorticity  at  that  time.  The  development  of  a  cyclostrophic 
balance  constraint  would  be  beneficial  for  flows  with  large  to  infinite  Rossby  numbers 
(Davis  and  Weisman  1994). 

Another  aspect  to  consider  follows  the  work  of  the  Gold  and  Nielsen-Gammon 
(2008  four  part  series).  The  used  PV  inversion  to  link  the  convective  parameters  to  the 
PV  distribution  based  on  the  invertibility  principle.  By  linking  features  like  the  cold  pool 
to  PV  anomalies,  the  effects  such  a  feature  has  on  the  convective  environment  might  be 
better  diagnosed. 
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While  not  related  to  the  avenue  of  this  research,  multiple  supercell  simulations 
were  conducted  with  various  microphysics  schemes.  A  large  difference  was  noted 
among  the  various  schemes  used  even  though  the  model  was  initialized  exactly  the  same 
each  time.  Since  each  of  the  microphysics  schemes  treats  hydrometeors  in  slightly 
different  way,  the  fonnation  of  certain  stonn  features  could  be  drastically  different.  This 
effect  could  be  seen  in  the  model  generation  of  low  level  rotation.  The  mid  level  storm 
rotation  acts  to  wrap  hydrometeors  (namely  rain)  around  the  updraft  while  the  storm- 
relative  wind  push  rain  away  from  the  updraft.  The  less  rain  available  to  wrap  around  the 
updraft,  less  evaporative  cooling  of  air  in  the  downdraft  would  produce  a  weaker  cold 
pool  and  weaken  the  low  level  baroclinic  generation  of  vorticity  along  the  cold  pool 
(Brooks  et  al.  1994a) 

More  diagnostics  are  also  need  to  verify  the  accuracy  of  the  inversion  technique, 
solving  for  the  omega  equation  to  compare  updraft  strengths  is  one  example.  Given  the 
good  qualitative  results  shown  here,  it  seems  feasible  that  PV  inversion  can  be  applied  to 
smaller  scales  than  what  it  was  initially  created  to  do. 
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